Beam Instabilities in the Scale Free Regime 
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The instabilities arising in a one-dimensional beam sustained by the diffusive photorefractive 
nonlinearity in out-of-equilibrium ferroelectrics are theoretically and numerically investigated. In 
the "scale-free model" , in striking contrast with the well-known spatial modulational instability, two 
different beam instabilities dominate: a defocusing and a fragmenting process. Both are independent 
of the beam power and are not associated to any specific periodic pattern. 



Introduction — Instability drives optical nonlinear waves 
into new and unexpected regimes. As for emergent be- 
haviors in complex systems, so instability can lead a sim- 
ple and featureless input wave to a whole spectrum of ef- 
fects that range from apparently random fragmentation 
to highly regular pattern formation. Systems support- 
ing solitons manifest a very general form of instability, 
modulational instability (MI), which fragments beams 
into patterns of periodic spots. [l|-|3] The entire process 
is rooted in the presence of one dominant spatial scale 
that forms in the nonlinear interaction, and this scale 
drives sufficiently wide waves into periodic patterns with 
a precise spatial period. 

Recently, the observation of a novel class of optical 
spatial solitons [H Q has been reported in disordered 
out-of-equilibrium ferroelectrics |l0L fll| . In this kind of 
media, the diffusion-driven photo-refractive nonlinearity 
12] can be largely enhanced to support "scale-free" non- 
diffractive beams, which can have arbitrary amplitude 
and waist. The mechanism is a direct by-product of the 
photorefractive band-transport model in disordered sys- 
tems, and applies to all photorefractive crystals, such as 
SBN and BaTiC>3, that are previously depoled into fer- 
roelectric clusters (see, for example, [13(), alongside or- 
ganic photorefractive glass [3] and photorefractive poly- 
mers [15(. More generally, scale- free optical effects stem 
from the interplay between diffraction and diffusion that 
mingle to cancel out their respective and intrinsic spatial 
scales, a condition that can occur in any wave support- 
ing systems where diffusion and diffraction both play a 
relevant role. One known example is in slow-light exper- 
iments, when atomic diffusion and paraxial diffraction 
accurately cancel each other out, depriving the wave phe- 
nomenon of a characteristic scale and cancelling diffrac- 
tion [3]. 

One-dimensional (stripe-shaped) scale-free solitary 
waves are observed to undergo a form of instability that 
breaks up (lCI |. a process reminiscent of MI. The puz- 
zling fact is that MI is wholly unexpected in a scale-free 
system, since, as we show below, the absence of a charac- 
teristic spatial scale provides no means for the propaga- 
tion to select a dominant perturbation. Remarkably, in 
agreement with experiments, we find that the scale-free 



regime, while curtailing the effects of MI, triggers entirely 
new beam instabilities. 

We theoretically investigate the dynamics of Gaussian 
scale-free one-dimensional solitons in two-dimensional 
propagation, and consider the spatial evolution of a pe- 
riodic perturbation with a transverse size different from 
that of the pump stripe beam. Although we find pertur- 
bations that can exponentially grow, these depend solely 
on pump beam waist and not on pump power, and no 
preferential spatial period emerges with maximum gain. 

We stress that the relevant one-dimensional nonlinear 
optical model (i.e., when neglecting one trasverse dimen- 
sion, as commonly done in MI theory), does not admit 
an unstable regime, and the considered beam breaking 
is the result of a subtle multidimensional coupling in the 
scale-free model, such that also the instabilities retain 
the scale-free character in the fact that no preferential 
spatial scale is developed during the dynamics, a feature 
potentially useful in imaging applications 17 1. 
Scale-free self-trapped beams — In the paraxial approxi- 
mation, for a linearly polarized beam with wavelength A, 
the slowly- varying optical field A (\A\ 2 — I is the opti- 
cal intensity) obeys the nonlinear propagation equation 
(hereafter denoted as the the scale-free model) [10 . Ill | 
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A = 0, (1) 



with k = 2nno/\, hq the bulk refractive index, and the 
characteristic length L given by 



L = 4imle y/gxPNR(K B T/q), 



(2) 



where Kb is the Boltzmann constant, T the tempera- 
ture, q the elementary charge, eo the vacuum permittiv- 
ity, Xpnr the low-frequency permittivity due to the polar 
nano-regions (PNR), and g > the relevant component 
of the quadratic electro-opitc tensor. 

As has been previously shown 1(J II |> Eq.([T]), admits 
bell-shaped soliton solutions as long as L > A. We specif- 
ically consider a perturbed striped beam. We begin writ- 
ing Eq.([T]) in dimensionless units as 
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where x = X/W , y = Y/W and z = Z/Z , with 
Zq = kW§ /2 the diffraction length, Wo an arbitrary 
beam waist, and u = L 2 /8X 2 , with a = 1/8 for L = A. In 
© "0 = A/ An is the dimensionless field scaled by an arbi- 
trary factor Ajv resulting from the intensity-independent 
feature of the scale-free model. 

The solitary waves of Eq.© represent diffraction- free 
beams with arbitrary intensity and waist [id . Ill| . Here 
we note that Eq. © also admits plane wave solutions for 
■ip = Aq = constant, which, at variance with Kerr me- 
dia, are not subject to any nonlinearly induced phase- 
shift (i.e., the nonlinear correction to the wave- vector is 
zero). We then consider the diffraction- free stripe (i.e., 
y— independent) solutions attained for a = 1/8 : 



ip s (x, z) = tpo(x) exp (i/3z) — Aq exp 



exp (if3z) . 



(4) 

with ipo(x) — exp(— x 2 /u> 2 ), (3 — —1/wq, and Aq and 
wq arbitrary independent parameters, such that this self- 
trapped beam exists at any intensity level and for any 
waist (within the validity of the paraxial approximation) . 
Note that as wo —> oo, /3 — > 0, one recovers the plane- 
wave solution ipo = Aq mentioned above. Solutions also 
exist for a > 1/8 (L > A) pUdl]], but nere we wm " focus 
on the case L = A (a — 1/8), as this is the condition 
achieved in experiments. 

Absence of Modulational Instability — We consider per- 
turbations to the exact solution, which is written as 



ip(x,y,z) = [tpo(x) + p(x,y,z)]exp(i/3z). 



(5) 



The linearized evolution equation for the perturbation p 
reads as 
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From Eq.©, one sees that if V>o = A (w a — > oo) no 
instability is expected (i.e., no solution such that p expo- 
nentially grows with z) ; the scale- free model of Eq. © is 
not exhibiting standard MI. This result is in striking con- 
trast with the well-known fact that MI always accompa- 
nies the existence of solitary waves solutions, as for Kerr, 
saturable or quadratic nonlinearities. Conversely, in 
the scale-free model, one has self-trapped bright beams, 
but no instability for the plane- wave solutions; unaffected 
by MI, the scale- free system turns out to be fertile ground 
for very different forms of instabilities, as we discuss be- 
low. 

Theory of the scale-free instabilities - We first note 
that Eq.© does not imply an exponential growth for a 
perturbation whose spatial profile in the x direction is 
the same as the pump beam ipo, i.e., for p such that 

p(x,y,z) = Mx)a+{zy kvV + Mx)u-(z)*e- lk * v . (7) 

Indeed the term containing p* in Eq.© disappears, so 
that there is no instability. Note that as Eq.© does not 



contain coefficients explicitly dependent on y, p can be 
expressed as a plane wave expansion with respect to y 
without loss of generality. We then write the perturba- 
tion as 

p(x,y,z) = 0i(a;)a + (z)e lfe « ! ' + ( /; 1 (a;)a_(z)*e- 4fc ^, (8) 

with ipi (x) an arbitrary spatial profile (specified below) , 
different from ipo(x). To keep the treatment as simple as 
possible, we limit analysis to the Gaussian soliton for L = 
X (a = 1/8), and we take the profile for the perturbation 
ipi (x) as a Gaussian with waist different from that of the 
pump ip (x), i.e., 
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such that [ipi,tj)i) = 1, with the scalar product (a,b) = 
J a*bdx. We use Eq.© in Eq.© and project over ipi, 
which corresponds to making an expansion in a Hermite- 
Gauss basis with respect to x and only retaining the first 
term of the expansion. We find the coupled equations for 
the amplitudes a± after Eq.©: 
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a± + ' u a T = 0. 

(10) 

Note that the last term coupling a± in Eq. ([T0l) is respon- 
sible for the instabilities and is proportional to w\ — w\ , 
hence for w\ = wq, we recover the result stated above, 
i.e., the absence of instability for a perturbation with the 
same a;— size of the pump beam. Analogously, the in- 
stability disappears for wo — > 00, corresponding to the 
plane- wave case, also discussed above. 
For wi 7^ wq, one finds that (fT0|) admits exponen- 
tially amplified solutions, which are written as a± = 
a± exp(Az), with the gain A given by 

4A 2 (fc v , m ) = {wl ~f or - ( e - X + ^ V ; (11) 



with r = wi/wq, one has 

4A 2 (fc„ r)w A = (r 2 - l) 2 - [(fc^ ) 2 + 1 - l/r 2 } 2 . (12) 

The most unstable perturbation corresponds to the val- 
ues r and k v that maximize A 2 . As detailed below, the 
analysis of Eq. (fT2"|) identifies two kinds of instabilities: 
with respect to perturbations of width greater than the 
pump beam (r > 1), and the opposite case (r < 1), 
denoted hereafter as defocusing and fragmenting insta- 
bilities, respectively. 

Defocusing instability — For a perturbation with w\ > 
wq, the condition A 2 > predicts the maximum gain at 
k y = and given by 



An = 



y/(r 2 - l) 2 - (1 - l/r 2 ) 2 



(13) 
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Ad is positive only for r > 1, and is shown in FigQJi. 
This maximum gain is not limited, and grows with r, 
thus revealing a self-propelling instability, such that if 
a perturbation with waist greater than the beam is su- 
perimposed, the beam tends to spread (the perturbation 
gains energy) and the spreading rate increases with the 
waist of the beam. Note that the gain is maximum at 
k y = 0, denoting an instability that does not tend to al- 
ter the striped shape of the beam by introducing period- 
ical modulations. This process is also more pronounced 
for small waists, as the maximum gain Ad goes like w^ 2 . 
We show in figure [2j an example of this instability, as 
obtained by numerically solving Eq.Q. The evolution 
reveals a defocusing of the beam, which is hence unsta- 
ble, and is compared with the linear case a = (linear 
propagation); for a > 0, the effect is more pronounced 
as the waist wo is reduced. Note that such an insta- 
bility is not observed for the two dimensional scale-free 
beams, neither in the numerical simulations nor in the 
experiments, [13] and these appear to be a sort of attrac- 
tor for the dynamics of the stripe beam, as shown in the 
following. 

Fragmenting instability — For r < 1 (perturbation 
smaller than the pump, i.e., w\ < Wo), the gain is maxi- 
mum at a k y > 0, fixed by r, and given by 
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The corresponding maximum growth rate A^ is 
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However, as r < 1, the maximum gain corresponds to 
r = (vanishing wi) with diverging k y , max , denoting the 
tendency of the beam to break up into very tiny spots, 
with no preferential spatial scale, in great contrast with 
the standard MI. Additionally, we note that, for a fixed 
r, the gain scales as the inverse squared waist, hence the 
more focused is the beam, the more pronounced is the 
instability. This is another remarkable difference with 
standard MI; e.g., in Kerr media as the pump power is 
increased, the gain grows as well; conversely in the scale- 
free model, the power does not affect the gain, which, on 
the contrary, increases when decreasing the beam spot 
size. Note also that the gain level for the fragmenting 
instability is lower than for the defocusing one and is 
limited by the upper value A_p(r = 0) = l/(2w 2 ), longer 
propagation distances are needed to appreciate its devel- 
opment. A notable outcome is that tiny details superim- 
posed onto the pump are amplified upon propagation. In 
Figj3^,b, we show an example of the fragmenting instabil- 
ity, by the evolution of a stripe perturbed by a Gaussian 
noise with 10% amplitude with respect to the pump. 

In Fig|3J:, we show the spectrum obtained numerically 
from Eq. , averaged over 10 noise realizations, and 
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FIG. 1: (a) Gain versus the ratio of the waist between the 
perturbation and the pump for r > 1 (maximum gain attained 
at k y — 0); (b) Left scale: gain versus r for r < 1; right scale: 
corresponding maximally amplified period. 
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FIG. 2: (Color online) Simulation of the defocusing insta- 
bility after Eq.([3|. (a) Output beam at z = 2 for a = 0.125, 
with a striped beam with wo = 1 and Ao = 1; (b) as in (a) 
with (7 = (linear propagation) ;(c) as in (a) with wo = 0.5; 
(d) beam waist in the x directions for (a), dotted line, for (b), 
dashed line, for (c) continuous line. 



compared with that expected from Eq. (|12[) : the shaded 
region corresponds to the overlap of the amplified spec- 
tral regions for various r and multiplied by a Gaussian 
spectral content. We report spectra for an initial waist 
wq = 0.5 and Wo = 4, showing that in the latter case 
the instability is moderated by about an order of magni- 
tude, as expected from Eq. (fT2)) . Spectra are arbitrarily 
shifted in the vertical axis for the sake of comparison. 
In recent experiments in copper doped KTN:Li, the de- 
scribed instability has been observed for with beam waist 
in the a;— direction of the order of 15/xm, after a propaga- 
tion distance of 6mm corresponding to Wo — 30/xm and 
w = 0.5 in FigH [Io| 

Conclusions — We have investigated wave instability 
in systems supporting scale-free optics, such as vitrified 
photorefractive dipolar glass. Without a dominant spa- 
tial scale, MI is found to disappear, outdone by two new 
forms of instability with a number of interesting phys- 
ical and mathematical properties. Scale-free instability 
is, in fact, able to amplify tiny beam perturbations at 
any spatial scale (in the paraxial approximation), with- 
out washing out or filtering the noise. The instability 
hence allows the wave to pick-up and carry large amounts 
of information, a fact that opens interesting perspectives 
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FIG. 3: (Color online) Simulation of the fragmenting insta- 
bility. The stripe solution is perturbed by a Gaussian noise 
p with 10% amplitude (\p/tpo\ = 0.1); (a) input beam; (b) 
output beam at z = 2 (Aq = l,wo = 0.5, a = 0.125); (c) 
spectrum obtained after the numerical simulations (continuos 
thick line for w$ = 0.5, thin line wq =4), the colored region 
is after eq. (|12|l showing the theoretically expected spectrum. 



for imaging applications. 
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